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Abstract. Graphene provides a fascinating testbed for new physics and exciting 
opportunities for future applications based on quantum phenomena. To understand 
the coherent flow of electrons through a graphene device, we employ a nanoscale probe 
that can access the relevant length scales - the tip of a liquid-He-cooled scanning 
probe microscope (SPM) capacitively couples to the graphene device below, creating 
a movable scatterer for electron waves. At sufficiently low temperatures and small size 
scales, the diffusive transport of electrons through graphene becomes coherent, leading 
to universal conductance fluctuations (UCF). By scanning the tip over a device, we map 
these conductance fluctuations vs. scatterer position. We find that the conductance is 
highly sensitive to the tip position, producing SG ~ e 2 /h fluctuations when the tip is 
displaced by a distance comparable to half the Fermi wavelength. These measurements 
are in good agreement with detailed quantum simulations of the imaging experiment, 
and demonstrate the value of a cooled SPM for probing coherent transport in graphene. 
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1. Introduction 

Graphene, a single atomic layer of carbon in a hexagonal lattice, has remarkable 
properties. It has conical conduction and valence bands that meet at a single point 
in fc-space (the Dirac point) [1]. Strong quantum confinement effects have been 
observed in quantum dots and nanoribbons |2], and the quantum Hall effect can be 
seen at room temperature |3J • Scanning tunneling microscopy has measured the surface 
topography local charge density and the local density of states (fl, 0|, and a 
scanned charge sensor has been used to map the ch arg e density 

Universal conductance fluctuations (UCF) occur when a coherent electron 



wave scatters repeatedly while it travels through a disordered conductor, following 
all possible paths through the sample. The different paths interfere with each other, 
creating a change in the conductance known as UCF that depends sensitively on the 
scatterer positions. When the size of the sample is less than the diffusive phase coherence 
length Ltf,, interference between paths yields a universal magnitude SG ~ e 2 /h for UCF, 



independent of the sample size and the degree of disorder. Theory [12|, |13J has predicted 
that the full UCF effect is obtained by moving a single scatterer a distance comparable 
to the Fermi wavelength A^. 

In this work, we use a liquid- He-cooled scanning probe microscope (SPM) (si. [l4 20 
to study coherent transport in graphene. We obtain conductance images that map the 
effect of a single scatterer on UCF. A charged SPM tip near the surface of a graphene 
sample creates an image charge that acts as a movable scatterer. This alters the electron 
wave function in the vicinity of the tip, leading to changes in quantum interference that 
give rise to UCF. An image of the sample conductance vs. tip position provides a spatial 
"fingerprint" that is unique to the arrangement of scatterers at a given Fermi energy. 
This technique allows us to observe the signatures of UCF without varying any external 
parameters (e.g. the magnetic field or gate voltage). Our approach reveals how UCF 
are created by the displacement of a single scatterer, as predicted by theory 
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UCF has recently been investigated in transport measurements of mesoscopic 



graphene samples |21l426j|. Our SPM technique provides a valuable spatial probe of 
coherent transport in graphene: 1) The tip can be adjacent to the two-dimensional 
electron gas, maximizing the spatial resolution, because graphene is two-dimensional 
material; 2) The Fermi energy Ep can be continuously varied from positive values for 
electrons, through zero to negative values for holes by using a back gate; 3) At T = 4 K, 
the observed coherence length (L^ ~ 500 nm) and elastic mean free path (l e ~ 50 nm) 
allow measurement in the coherent regime. Our results for graphene should also apply to 
other two-dimensional conductors, though some questions may be raised by scattering 
in graphene's unusual band structure 
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Figure 1. (a) Schematic of the graphene sample mounted in the cooled scanning probe 
microscope (SPM), showing a Hall bar contacted by six Cr/Au leads. The SPM tip and 
lead 4 are grounded, with a 25 nA rms current between leads 1 and 4 at 5 kHz. Voltage 
is measured between leads 2 and 3 using a lock-in amplifier. A back-gate voltage V g is 
applied to the degenerately doped Si substrate, (b) Measured conductance G vs. V g , 
with a linear fit for V g < 20 V. (c) Three consecutive measurements of G vs. V g with 
the tip fixed far from the sample, with the linear background from (b) subtracted, (d) 
Same as (c), but with the tip 10 nm above the sample, at three different locations 
spaced 100 nm apart. 



2. Experimental Methods 

The graphene samples studied in these experiments are single-atomic-layer Hall bars, 
with a geometry shown schematically in figure UK The experimental data in the figures 
below are from a sample with width 500 nm and voltage contacts (leads 2 and 3 in 
figure UK) with centers spaced 1200 nm apart. Graphene flakes were prepared through 
mechanical exfoliation (the "sticky tape method") and deposited onto a degenerately 
doped Si substrate capped with 280 nm of Si0 2 . A back gate voltage V g is applied 
between the substrate and the graphene to vary Ep and the charge density n. Using 
electron beam lithography, Cr/Au leads are deposited onto the graphene, after which 
the Hall bar structure is formed via an oxygen plasma etch. The presence of single- 
layer graphene is confirmed by observing quantum Hall conductance plateaux at the 
expected values 4(z/ + l/2)e 2 /h, where v is an integer. The sample is then mounted 



on a home-built scanning gate microscope |15l . Il9j and cooled in He exchange gas in 



thermal contact with a liquid He bath at T = 4.2 K. Similar results have been obtained 
on three other samples with similar dimensions. The four-probe conductance G, shown 
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in figure [lb, displays the characteristic linear variation of G vs. V g on either side of the 
Dirac point at Vbwoc — 22 V. 

Coherent, diffusive transport is expected when the sample size L < L^, the 
electron's diffusive coherence length, and L 3> l e , the elastic mean free path. At 
T = 4 K, we obtain « 0.5 /mi from weak localization and magnetoconductance 
measurements (not shown). is larger than the 0.4 x 0.4 /im 2 field of view in figure |2l 
and comparable to the sample width (W = 0.5 /im) and length (L = 1.2 /im); all of 
these lengths are much larger than l e ~ 50 nm. 

From the slope of G vs. V g in figure [lb and the capacitance between the back gate 
and the graphene, discussed below, the electron and hole mobility away from the Dirac 
point is found to be fi ~ 7200 cm 2 /Vs. The shift of the Dirac point from V g = to 
Voirac = 22 V is attributed to charged impurities located either above or below the 
graphene layer, which induce a charge in the graphene. 

To create a movable scatterer, a conducting, voltage-biased SPM tip with radius 
of curvature r t i V = 20 nm is held at a height hu p = 10 nm above the graphene. In the 
measurements presented here, the tip is grounded, so that the charge on the tip is set 
by the contact potential between the tip and the graphene. In addition, image charges 
are created in the tip from impurities on the surface of the graphene sample. 

The spatial profile of the density perturbation created in the graphene by the 
SPM tip was computed using classical electrostatic finite-element simulations (Maxwell, 
Ansoft LLC). The graphene is modeled as a planar conductor, with the observed 
offset Voirac of the Dirac point modeled by a homogeneous layer of charge above the 
graphene. The tip is realistically shaped and located above the sample at a height 
hti P = 10 nm. The back gate is modeled as an infinite conducting plane, separated from 
the sample by 280 nm of Si02- The average carrier density in the graphene is found to 
be n = a(V g — Voirac), with a = 8.3 x 10 10 cm -2 V -1 . The spatial profile of the image 
charge created by the SPM tip in the graphene layer has a maximum ~ 5 x 10 11 cm -2 
and a Lorentzian-like shape with half-width at half maximum (HWHM) r scat « 25 nm. 
The size and magnitude of the tip perturbation can be compared to the naturally 
occurring variations in charge density (charge puddles) in graphene, which are found 
experimentally Sag and theoretically Qto have charge densities ~ 4 x 10 cm 
with a characteristic diameter ~ 20 nm, for graphene flakes on a Si02 substrate. The 
perturbation to the charge density created by the SPM tip has approximately the same 
amplitude, and about double the spatial size of these pre-existing inhomogeneities. 

The conductance fluctuations visible in figure [Tb can be identified as UCF. They are 
reproducible, and have a root-mean-squared (rms) magnitude 5G = 0.64 e 2 /h. FigureHfc 
shows the conductance fluctuation AG vs. V g when the tip is fixed far from the sample 
with tip height h tip > 100 /im. The three traces from consecutive V g sweeps show good 
reproducibility. Bringing the charged tip near the graphene ihu v = 10 nm) creates an 
image charge in the electron gas that significantly alters the conductance fluctuations. 
Three AG vs. V g traces in figure [TH for different tip positions spaced 100 nm apart, 
demonstrate that UCF is sensitive to the spatial configuration of scatterers - a change in 
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Figure 2. (a) Conductance image G(r) vs. tip position r at T = 4 K, for tip 
height h Up = 10 nm. The density n = -2.7 x 10 12 cm" 2 (V g = -10 V) is far 
from the Dirac point. The 400 x 400 nm 2 scan area is located in the center of 
the sample, as indicated in the schematic diagram, (b) Same as (a) except for a 
density n = —1.2 x 10 11 cm~ 2 (V g — 20 V) near the Dirac point, (c) Repeatability 
is demonstrated by two 400 x 400 nm 2 conductance images taken in succession, a few 
minutes apart, with the rightmost panel showing the difference between the first two 
images. 



the position of a single scatterer is enough to decorrelate the conductance fluctuations, 
as predicted by theory [12, ll3 . 



3. Results 



3.1. Experimental UCF images 



Using our SPM, we can study UCF by controllably raster scanning the tip position 
over an area of the sample; previous studies of the effect of single scatterers were 

Figures [2^ and [2b show 

and 
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based on charge hopping at random positions 

conductance images of G vs. tip position r at densities n = —2.7 x 10 



12 cm 2 



n = —1.2 x 10 11 cm -2 respectively, in a 400 x 400 nm 2 area located at the center of 
the sample. At high density (figure EK), conductance fluctuations are observed with 
rms magnitude SG = 0.35e 2 /h and characteristic lateral size ~ 10s of nm in agreement 
with UCF theory, as shown below. At low density, near the Dirac point (figure Eb), the 
conductance fluctuations have a smaller magnitude 5G ~ 0.1 e 2 /h and a larger lateral 
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Figure 3. (a) Conductance images G(r) vs. SPM tip position r in a 400 x 400 nm 2 
region of a graphene sample at different back gate voltages V g . The color scale spans 
a range of ± le 2 /h. An arrow points to the same location in each image, highlighting 
their continuous evolution, (b) A repetition of the series of images in (a), performed 
1.5 hrs later, demonstrating the repeatability. The arrows point to the same location 
and feature as in part (a). 



size ~ 100 nm. Previous transport measurements show the magnitude of UCF in single 
and multilayer graphene 2l|, |24], |26( decreases monotonically towards the Dirac point. 
The conductance images shown in figure 2 are reproducible, as expected for UCF. Two 
images taken ~ 3 min. apart (figure [2b) are nearly identical; the difference shows only 
small changes, likely caused by the motion of charged defects in the substrate. The 
UCF images shown in figure [2] are loosely analogous to speckle patterns produced by 
the coherent scattering of light in a diffusive medium. 

Another test for UCF is obtained by measuring the correlation between conductance 
images recorded at different Fermi energies Ep by changing the density n with the back 
gate voltage V g . The correlation energy E c is the range of Ep over which UCF remain 
correlated. We can determine E c from our UCF conductance images by finding the 
change in V g needed to reduce the correlation between two images by one half; this is 
the correlation voltage V c . 

Figure [3b, shows a series of conductance images recorded at backgate voltages 
decreasing from V g = 0.8 to 0.2 V in 0.1 V steps. By eye, one can see that the images 
evolve smoothly from one to the next, becoming less correlated as the change in V g is 
increased. For example, the arrows point to the location of a dark spot at V g = 0.5 V in 
the middle of the series, which has almost completely disappeared at V g = 0.2 or 0.8 V. 
The reproducibility of these UCF images over time intervals ~ 1 hr is demonstrated by 
figure [3b which shows a repetition of the V g series, performed 1.5 hrs later. The arrows 
point to the same location as in figure EK, showing that the same feature remains. 

The correlation Cab between two conductance images Ga (r) and G7,e(r) vs. SPM 
tip position r, is Cab — J (GU( r ) ~~ (GU))(G\B(r) — (G7 B ))<ir, where angle brackets denote 
the average over r. The normalized correlation Cab, such that the autocorrelation of 
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an image is equal to unity is 

Cab 



C 



AB 



From a series of conductance images at different back gate voltages V g , we obtain 
the normalized correlation C(v g )(v g +AV) between two images, GV 9 (r) and GV 3 +Av( r )) 
separated by a fixed change AV in V g . The average correlation (C(v g )(v g +AV))v g vs - AV 
is then obtained by averaging over different values of V g for a fixed AV. 
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Figure 4. (a) Normalized correlation (C(y )(y +Av))v g vs - between two UCF 

conductance images recorded at different back gate voltages V g and V g + AV, averaged 
over V g . The correlation voltage V c is the half-width-at-half-maximum (HWHM) of 
this curve, (b) Correlation voltage V c vs. V g . Points: V c computed from experimental 
conductance images. Line: theoretical curve (see text) following Ref. [jjj ). 



Figure UK shows the normalized correlation (C(y g 



)(V g +AV)/V g 



)v g , averaged over V 6 



9 1 



VS. 



AV between two UCF images recorded at different back gate voltages V g and V g + AV 
from V g = —1.0 to 1.0 V in steps AV = 0.1 V; the UCF images are similar to those 
in figure |3K As shown, the correlation falls off as AV moves away from zero. The 
correlation voltage V c is defined as the halfwidth-at-half-maximum (HWHM) of this 
curve. To find 14 at a particular back gate voltage V g = V g , a series of conductance 
images is recorded over the range V„ = V° — 1.0 V to V„ = V° + 1.0 V in steps of 0.1 V. 
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Figure 5. (a) Simulated UCF conductance image G(r) vs. tip position r in graphene 
for a density n = —8 x 10 11 cm~ 2 away from the Dirac point, and (b) simulated UCF 
conductance image for a density n = —4 x 10 10 cm~ 2 near the Dirac point. These 
simulations show a magnitude SG and spatial size of UCF comparable to the data in 
figure [2j The difference in conductance G between the simulation and experiment is 
caused by the difference in sample aspect ratio and mobility. 



We then calculate (C(y g )(v g +Av))v B vs - for this series of images, and the correlation 
voltage V c is obtained from the HWHM of this curve. For the data in figure we find 
V c = 0.27 V at V g = 0. Note that the range of V g (2V) covered by this procedure is 
small compared to the full range of V g (10s of V) considered in the experiments. 

A plot of the measured correlation voltage V c vs. back gate voltage V g 
is shown by the red dots in figure HJo. Theoretically, the correlation voltage 
V c = (2E c /hvo)(\V g — VDiracl/Tta) 1 / 2 where E c is the correlation energy, and vq = 
1.1 x 10 6 m/s is the Fermi velocity. The correlation energy is E c 2.8 ksT in the 
thermally broadened regime 25] that is appropriate here. The theoretical curve for V c 
at T = 4 K, is shown by the solid curve in figure Hb; it is in good agreement with 
the experimental results providing strong evidence that our conductance images display 
UCF. 

To ensure that the measured loss of correlation with increasing change AV in back 
gate voltage V g is a repeatable effect and not caused by random drift over the course 
of the measurement, we calculate the correlation Cman between a conductance image 
Gv g ,t{ r ) obtained at time t and the same scan Gy t /(r) repeated at t' — t + 1.5 hrs, 
as shown in figures [3|i and [3b- The correlation for different times (C(t)(t'))v g ~ 0.5, 
averaged over V g , is much higher than the correlation (C(v g )(v g +AV))v g f° r different back 
gate voltages separated by A V = ±1 V, as shown in figure Hk, demonstrating that the 
complete loss of correlation with gate voltage cannot be caused by drift. 
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3.2. Simulated UCF images 

We have performed quantum simulations of coherent transport in graphene including 
the potential from a movable tip (TsJ , shown in figure EJ to clarify the origin of features 
seen in the measured conductance images. The simulation results display UCF that 
change with the tip position, producing conductance images in good agreement with 
the experiment. The simulated UCF images show the same trends with Ep as in 
the experiment. By varying parameters in the simulation, such as the size of the tip 
perturbation, we can further explore the effect of a movable scatterer in graphene, and 
compare to the predictions of analytic theories. 

Theoretical simulations of UCF conductance images for graphene were obtained 
by using a finite-difference method to calculate the conductance through a disordered 



graphene sample [31[ . The potential in the graphene layer is given by the combination of 
a local potential from the tip and a number of randomly placed electrostatic scatterers. 
Details of the simulations are given in the Appendix. 

By calculating the conductance G(r) vs. the tip position r, we produce the same 
type of conductance images as obtained in the experiments. Theoretical conductance 
images, such as figure El are obtained by rastering the tip position r in a plane above 
the sample. The conductance G at each tip position r is simulated using the combined 
potential U from the tip, the intrinsic scatterers and the back gate. Each image consists 
of 80 x 80 evaluations of G, with tip positions spaced 5 nm apart, centered within the 
sample area. 

The simulated images shown in figure[5]are in good agreement with the experimental 
results in figures [2^ and [2b. Figure EK shows a simulated conductance image at a density 
n = —8 x 10 11 cm -2 far from the Dirac point. The image displays spatial conductance 
fluctuations with amplitude 5G ~ e 2 /h, and lateral size ~ 10s of nm, similar to the 
fluctuations observed in the experiment in figure Close to the Dirac point (figure Eb, 
n = —4 x 10 10 cm -2 ), the simulated images show UCF with larger lateral size, in 
agreement with figure [2b. The simulated UCF images have characteristics similar to 
the measured images: rms amplitude SG ~ e 2 /h, and spatial size of UCF that increases 
near the Dirac point. The agreement between the simulations and the experiment verify 
that the measured images show UCF caused by the motion of a single scatterer. 



3.3. Spatially-resolved UCF 

Our measurements provide a unique ability to probe theoretical predictions 0, 
for the effect of a single movable scatterer on UCF. To quantify the spatial size of 
the features in a conductance image G(r), we calculate the spatial autocorrelation 
C(tq) = J G(r)£r(r— ro)dr. FiguresEh andEb show C(tq) away from the Dirac point for 
the experimental and simulated results, respectively. The width of the central peak in 
these plots corresponds to the spatial size of the fluctuations in the original conductance 
image. Figures Eb and EH show C(r ) for Ep close to the Dirac point, with significantly 
broader peaks. 
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Figure 6. (a) Experimental (exp) and (b) simulated (sim) autocorrelation C(ro) of 
UCF conductance images away from the Dirac point [(a) n = —7.1 x 10 11 cm -2 and 
(b) n = —8 x 10 11 cm -2 ]. White = high and black = low. (c) Experimental (exp) 
and (d) simulated (sim) C(r ) near the Dirac point [(c) n = — 8 x 10 10 cm~ 2 and (d) 
n = —4 x 10 10 cm' 2 ]. Diagrams of graphene band structure schematically indicate 
the Fermi energy Ep in (a) to (d). (e) Simulated correlation length l c at a fixed 
n = —8 x 10 11 cm~ 2 vs. the radius r scat of the movable scatterer created by the tip, 
with linear fit (red). Points are the average of 5 disorder configurations. Dotted line 
shows the analytical prediction l c — 0A6Xp. (f) Measured correlation length l c vs. V g 
(points) obtained from autocorrelations C(tq), shown with the analytical prediction 
l c = 0.46Af (dotted line), and an empirical fit (red line) to l c — 0.46Ai? + ro with 
ro = 22 nm; the value of ro obtained from the fit is close to the SPM tip radius, as 
expected. Data points represent the average of four images at slightly different V g to 
reduce noise. 



We find that the correlation length l c for a UCF image is determined by the Fermi 
wavelength. The correlation length can be extracted from C(yq) by averaging over the 
angular dependence and defining l c to be the HWHM of the resulting curve. Figure |6]F 
shows l c vs. V g from a series of experimental conductance images spanning the Dirac 



point. Theory 12] predicts l c ~ 0.46Ai? (dotted line in figure [6f), where the Fermi 
wavelength A^ = 2(7r/|n|) ' 5 . Both theory and the data show a peak in l c at the Dirac 
point. The analytical calculation, however, does not take into account the spatial extent 
of the tip or the long-range scatterers (10s of nm) which are not negligible. 

The effect of the spatial size r scat of the movable scatterer created by the tip is 
investigated in simulated images by finding l c vs. r scat , shown in figure |6fe, where 
r scat is the HWHM of the image charge density. The scatterer effectively smears 
out the fluctuations on a length scale (1.4 ± 0.3) x r scat , shown by the red line in 
figure Ek. This smearing can be modeled by adding an offset ro to the correlation length, 
l c = 0.46Ai? + ro; a best fit shown by the red line in figure Ef is ro = 22 ± 1 nm, which 
corresponds closely with the tip radius r tip w 20 nm and with the size r scat w 25 nm of 
the image charge created by the tip. 
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4. Conclusions 

Our SPM imaging technique probes how coherent transport through a mesoscopic 
graphene sample is affected by the motion of a single scatterer. By scanning a charged 
SPM tip over a graphene device, we obtain conductance images that display fluctuations 
with amplitude 5G ~ e 2 /h and spatial size ~ 10s of nm comparable to the Fermi 
wavelength. The UCF conductance images repeat over times up to 1.5 hrs, as predicted 
for UCF. The correlation between two images is destroyed by changing the Fermi energy 
and density, by changing the back gate voltage V g . The correlation voltage V c obtained 
from our measurements agrees well with the theoretical prediction for the correlation 
energy for UCF. It is striking that a change AV g < 1 V in back gate voltage is sufficient 
to completely change the conductance images; this is expected for UCF created by the 
interference of electron waves traveling along different paths, and demonstrates that the 
images are not simply reflecting the underlying charge density puddles. We see that the 
interference that gives rise to UCF is highly sensitive to the position of even a single 
scatterer, yielding the full fluctuation 5G ~ e 2 /h when the tip is displaced by only 
several 10s of nm. 

To verify that the observed conductance images represent UCF caused by the 
motion of a single scatterer, we simulated the effect of the tip-created scatterer on 
the conductance of a graphene sample. The simulated conductance images reproduce 
the features seen in the experimental images: fluctuations 5G ~ e 2 /h with lateral size 
~ 10s of nm, which depend sensitively on the Fermi energy and the arrangement of the 
scatterers. The simulations also confirm the observed increase in the spatial size of the 
fluctuations near the Dirac point. 

Because universal conductance fluctuations result from quantum interference, one 
would expect their spatial length scale to depend on the electron wavelength. Indeed, 
from the experimental conductance images and numerical simulations, we find that the 
spatial size of the fluctuations is comparable to the Fermi wavelength Xp. We obtain 
good agreement with theoretical predictions, taking into account the realistic spatial 
size of the tip-created scatterer. 

Our measurements demonstrate the utility of a low-temperature scanning probe 
microscope for studying the coherent flow of electrons through graphene. The 
conductance images shown above provide a spatial view of how the interference of 
electron waves leads to UCF. This imaging technique will also be useful for the 
investigation of magnetoconductance fluctuations and weak localization in graphene. 
By using a probe of size comparable to the electron wavelength, we gain new insight 
into the quantum behavior of electrons as they flow through a graphene device. 
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Appendix: Simulation methods 

The numerical calculations we have performed to model our results follow the method 



described in Ref. [31j. The sample is discretized into a square lattice and the Dirac 
equation is solved using a finite difference method on this grid. An ideal lead is connected 
to both sides of the sample with propagating electron modes incident on the sample 
edges. The conductance is obtained by calculating the transfer matrix for these modes 
as they travel across the sample. Further details of these simulations will be given in a 
separate publication. 

In our simulations, the sample grid consists of 102 x 153 points, spaced 5 nm apart. 
This places a lower limit on the Fermi wavelength of \ m in = 10 nm, corresponding to a 
maximum charge density of n max = 47r/A 2 rein = 4n x 10 12 cm -2 . The direction of current 
flow is across the narrow dimension of the grid (510 nm), with a width of 765 nm in the 
transverse direction. Periodic boundary conditions are applied at the transverse edges, 
and we focus only on the 400 x 400 nm 2 square in the middle to avoid effects of the 
boundary conditions. Note that the aspect ratio of the simulated sample L/W = 2/3 
is less than the aspect ratio L/W = 2.4 for the experimental sample. For the same 
conductivity, the conductance G for the simulation will be a factor ~ 3.6 larger than G 
in the experiments. 

Disorder in the graphene is modeled in the simulations as a sum of screened 
electrostatic potentials created by point charges located above or below the graphene 
layer. According to the method of images, a point charge q located a height a above a 
conducting sheet induces a charge density in the sheet 

= 27 r(p 2 Ta 2 )3/ 2 (A ' 1} 

where p is the radial coordinate away from the position of the point charge. We 
then build up the total disorder charge density cr^x) vs. position x as a sum of 
such functions, centered at randomly chosen lattice sites, with a fraction n« = 0.2 
of lattice sites occupied. The charge q = 2.5 e is chosen to yield a rms charge 
density ~ 4 x 10 11 e/cm 2 that is in agreement with the observed charge puddles 
in scanning tunneling j^, 0] and scanning charge sensor measurements The sign of 
each impurity charge is randomly chosen to be positive or negative, with equal numbers 
of positive and negative charges. (The offset of the Dirac point from V g = in the 
experiment is not explicitly included in the simulation.) The distance of the charged 
impurities from the sample is set to a = 10 nm to match the lateral size scale of 
the image charge puddles with puddles observed in scanning tunneling and scanning 
charge sensor experiments BBS. These numbers combine to yield an effective density 
n<imp = 2 x 10 12 cm -2 of impurities with charge ±e. The resulting simulated conductance 
increases linearly with n away from the Dirac point, and is rounded off at n = 0. The 
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mobility of the simulated sample is // ~ 15000 cm 2 /V s, which is a factor of ~ 2 
greater than the experimental mobility. The minimum conductivity is 5.3 e 2 /h in the 
simulations, in agreement with the measured value 5.7 e 2 /h. 

We simulate a conductance image by adding to the charge density an additional 
perturbation a t i P created by the tip, centered at position r. We model the tip as a point 
charge above the sample, so a tip has the functional form given in equation IA.1I We 
adjust the tip height a to control the width of the tip perturbation r scat . For the images 
in figure El we chose a = 10 nm, the same distance from the sample as the charged 
impurities. To test the effect of the radius of the scatterer created by the tip, we vary 
a = 10 nm to a = 32.5 nm in figure [6fe. We set the tip charge q to yield a peak image 
charge density a max ~ 5 x 10 11 e/cm 2 , as determined from electrostatic simulations, 
described above. 

Finally, an overall offset o"o to the charge density is added to yield the desired Fermi 
energy, controlled by V g in the experiment. The total charge density in the graphene is 
then given by a = a + ad + o- tip . Using the relationship between the Fermi energy and 
charge density in graphene, we can now find the potential vs. position x in the graphene 
layer: 

U(x) = hv x sgn(o-(x)) v / 7r|cr(x)| (A.2) 

where vq is the Fermi velocity, and "sgn" is the sign function. This potential is 
then plugged into the simulation to model the disordered potential through which the 
electrons flow. Note that the square root in equation IA.2I means that the different 
contributions to the potential do not add arithmetically to the total potential. That is, 
the contribution to the potential from disorder becomes smaller as the overall charge 
density increases. 

The empirical model for the disorder described above is based on calculations 



for charge puddles caused by screened impurities in graphene in Refs. [28|, |32|, and 



measurements of the size and magnitude of charge puddles in Refs. jsl, Q, 0] . We neglect 
short-range scattering with lattice defects, scattering from ripples or trig onal warping, 
and the quantum corrections to the screening expected at low density [32|. We find that 
this simple model of long-range, ideally-screened electrostatic scatterers is sufficient to 
reproduce and understand the experimentally observed phenomena. 
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